Frequency evaluation of the doubly forbidden 1 So — > 3 Po transition in bosonic 174 Yb 
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We report an uncertainty evaluation of an optical lattice clock based on the So ^ 3 Po transition 
in the bosonic isotope 174 Yb by use of magnetically induced spectroscopy. The absolute frequency 
of the 1 So <-> 3 -Po transition has been determined through comparisons with optical and microwave 
standards at NIST. The weighted mean of the evaluations is i/( 174 Yb)=518 294 025 309 217.8(0.9) Hz. 
The uncertainty due to systematic effects has been reduced to less than 0.8 Hz, which represents 
1.5 x 10~ 15 in fractional frequency. 

PACS numbers: 32.10.Dk, 06.30.Ft, 32.70.Jz, 39.30.+W 



Optical frequency references are now reaching frac- 
tional frequency uncertainties lower than those of the 
best Cs primary standards, which define the second. In- 
tense efforts to reduce these levels to well below one part 
in 10 16 are underway with many different atomic sys- 
tems (neutral and ionic) for a variety of applications in- 
cluding tests of fundamental physics and the develop- 
ment of the next generation of primary frequency stan- 
dards [H, 0]. Here we report the first detailed evalua- 
tion and absolute frequency measurement (with an un- 
certainty of 1.5 x 10~ 15 ) of a clock system based on neu- 
tral 174 Yb atoms confined to a one-dimensional optical 
lattice. Ytterbium is an attractive atom for optical clock 
studies due to its large mass and numerous abundant iso- 
topes that offer a variety of nuclear spins (0,1/2,5/2). For 
this work, we focus on the potential of spin-zero isotopes 
for lattice clock work with the first evaluation of such a 
system at the 1 x 1CP 15 level. 

Spin-zero isotopes are appealing for frequency metrol- 
ogy due to their simple structure. Indeed atomic clocks 
based on J = <-> J' = transition in these iso- 
topes present a nearly pure two-level system in that it 
lacks Zeeman structure, first-order sensitivity to mag- 
netic fields, optical pumping effects, and first-order vec- 
tor/tensor sensitivity to the lattice laser intensity and 
polarization. Moreover, the spin-zero (even) isotopes in 
the alkaline earth-like atom systems are generally more 
abundant and easier to prepare (through laser cooling 
and trapping) for high resolution spectroscopy. To realize 
these advantages, however, it is necessary to use an ex- 
ternal field such as a magnetic bias field (as in this work) 
to induce a nonzero transition strength for the otherwise 
truly forbidden transition [H, 0| ■ The presence of this ex- 
ternal field leads to a small shift of the clock frequency 
that must be calibrated. Additionally, the weakness of 
the induced transition requires higher probe light inten- 
sities than are used with the odd (half-integer spin) iso- 
topes, which lead to an AC Stark shift that also requires 
careful evaluation. This is the essential trade-off between 
the odd and even isotopes: one accepts two larger system- 



atic effects in exchange for a considerably simpler atomic 
system that does not require optical pumping or hopping 
between multiple spectroscopic features in order to sup- 
press first-order magnetic field and lattice polarization 
sensitivities. Optical lattice clock studies with Sr have 
thus far been focused on the odd isotope ( 87 Sr, with nu- 
clear spin 9/2), and the one evaluation with an even iso- 
tope ( 88 Sr) had an uncertainty considerably higher than 
those of its odd counterparts [5j . Here we show that the 
even isotopes, especially in the case of Yb, which has 
somewhat more favorable atomic properties for even iso- 
tope work [3[, can achieve extremely low uncertainties 
as well, with no identified barriers to reaching the 10~ 17 
level. 

Experimental details on the cooling and trapping setup 
have been reported in Q. In brief, 174 Yb atoms are 
trapped and cooled in a two-stage magneto-optical trap 
(MOT) working on the strong dipole allowed 1 So <^> 1 P 1 
transition at 398.9 nm and then on the intercombina- 
tion 1 S'o <-> 3 Pi transition at 555.8 nm. Typically 10 5 
atoms are trapped at a temperature of less than 40 fiK in 
350 ms. About 10 4 atoms are then loaded to a ID stand- 
ing wave optical lattice operating at the magic wave- 
length of 759.35 nm for the 1 So ^ 3 Pq clock transition. 
A lattice potential depth of about 500 E r (recoil energy 
Er/ks = 100 nK for Yb at the lattice frequency) is real- 
ized by tightly focusing (waist « 30 /im) a 1 W injection- 
locked ThSapphire laser. After a 25 ms delay that allows 
for the quadrupole magnetic MOT field to dissipate and 
the static magnetic field to be turned on for the spec- 
troscopy, the 1 S'o <-> 3 Po clock transition at 578 nm is 
probed with a 7r-pulse lasting 50 to 150 ms. The probe 
laser is collinear with the lattice beam and is focused 
onto the atoms with a waist of 60 /im. Both the lattice 
polarization and the probe laser polarization are aligned 
with the vertical static magnetic field. 

The 578 nm clock laser is based on sum-frequency 
generation of a Nd:YAG laser at 1.32 /im, and an Yb 
fiber laser at 1.03 /im in a periodically-poled lithium 
niobate waveguide [7|. The 578 nm light produced is 
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then frequency stabilized with a Pound-Drever-Hall lock 
to a vertically-mounted high finesse stable cavity. The 
short-term frequency stability is 2.2 x 10~ 15 at 1 to 2 
s, with a residual frequency drift of less than 0.4 Hz/s. 
The residual drift is monitored and canceled to first or- 
der through feed-forwarding to an acousto-optical mod- 
ulator interposed between the laser source and the ref- 
erence cavity. The stabilized yellow light is then sent 
through phase-noise-compensated fibers to both a fem- 
tosecond frequency comb for clock comparison and to 
the atom trap. For a typical magnetic field strength of 
B = 1.3 mT and a probe intensity of Iq — 300 mW/cm 2 , 
the Rabi frequency is Q ~ 5 Hz, which generates a 
Fourier-limited linewidth of about Av = 10 Hz. To pre- 
vent spurious ac Stark shifts, the pre-cooling and trap- 
ping beams are switched off during the spectroscopy with 
acousto-optical modulators and mechanical shutters. 

The amount of population transfer via the clock transi- 
tion is determined through fluorescence detection on the 
strong 399 nm transition of the atoms remaining in the 
1 5'o ground state. To lock the frequency of yellow probe 
laser to the atomic transition, an error signal is derived 
by alternately probing the half-width points and then 
demodulating with a digital microprocessor, which pro- 
vides frequency corrections to an acousto-optical modu- 
lator. For a typical observed S/N ratio of ~ 10 and a 
probe duty cycle of 15 %, the projected clock stability is 
1 x lO-^r-V 2 . 

Evaluation of the systematic shifts of the Yb lattice 
clock was accomplished through comparisons with other 
optical frequency standards at NIST. While initial mea- 
surements were performed against a neutral Ca stan- 
dard , the measurements reported here were performed 
against the Hg + optical standard [9j. For the compar- 
isons an octave-spanning TkSapphire optical frequency 
comb 10] was locked to the Hg + standard, and the beat 
frequency between the comb and the Yb standard was 
counted with 1 s gate times. We typically changed the 
value of the experimental parameter under study every 
200 to 300 s in order to determine the eventual shift with 
a statistical uncertainty of about 3 x 10~ 16 (see Fig. [5]). 

In the case of spectroscopy of the clock transition 
with the even isotope, the linear Zeeman effect is zero, 
with only a second-order dependence on the magnetic 
field amplitude. In Fig. []Ji we show an evaluation of 
the second-order Zeeman shift obtained by varying the 
static magnetic field between 1 mT and 3.5 mT. The 
resulting second-order Zeeman coefficient in 174 Yb is 
(3 = —6.12(10) Hz/mT 2 , found by a quadratic fit to zero 
field. For a field value of 1.6 mT the shift is approxi- 
mately 18 Hz and is determined with an uncertainty of 
less than 0.3 Hz. The uncertainty in the systematic shift 
(Ab = 2/3|B|i5B) is limited mainly by the knowledge of 
the applied static magnetic field B, which is calibrated to 
better than 10 /xT through spectroscopy of the magnet- 
ically sensitive Sq <-> 3 Pi transition. Strategies for re- 
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FIG. 1: Systematic shifts of the 1 S'o «-» 3 fb transition in 
Yb optical lattice clock, a) Measurement of the second or- 
der Zeeman shift. The second order coefficient is evaluated 
/3 = -6.12(10) Hz/mT 2 . b) Probe ac linear Stark shift. The 
coefficient k = 15(3) mHz/(mW/cm 2 ) is in agreement with 
the theoretical estimation (see text), c) Density shift. The 
data points are consistent with zero density shift with a to- 
tal deviation of 0.6 Hz. d) Servo error. The frequency offset 
due to uncompensated linear drift of the reference cavity is 
shown (slope = 1.82(10) Hz/(Hz/s)). With periodic manual 
correction of the feed-forward system used to remove the lin- 
ear cavity drift this uncertainty can be made negligible. 



ducing the uncertainty of the second-order Zeeman shift 
include refining the absolute calibration of the magnetic 
field down to less than 1 //T (one possibility could be the 
use of the narrower Sq ^-> 3 -Fb transition in the fermionic 
171 Yb isotope) and operation at lower bias field values. 

In a similar way, the ac Stark shift induced by the yel- 
low probe light has been determined by measuring the 
shift as a function of intensity and extrapolating to zero 
intensity. Although the probe laser power is controlled 
with good precision (< 1%) with a photodiode and fast 
servo control, the absolute intensity at the location of 
the atoms is not well determined. In the current setup, 
the clock excitation light is focused to a ~ 60 /im waist 
at the atoms, which is larger by only about a factor of 
two than the waist size of the lattice beam. This im- 
plies that the atoms experience an average intensity and 
not uniform illumination. In addition, drifts in alignment 
can cause the average intensity experienced by the atoms 
to vary by 10 % from day to day in the current setup. 
Fortunately, the intensity experienced by the atoms is 
directly proportional to the optical power, so the shift 
due to the yellow probe can be accurately extrapolated 
to zero intensity. As is reported in Fig. |T|b, this effect 
shows a linear dependence on the probe intensity with 
a slope of k ~ 15(3) mHz/(mW/cm 2 ), in agreement 
with the theoretical estimate [3J. In this case, the typi- 
cal ^300 mW/cm 2 probe intensity gives a total shift of 
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about 7 Hz with an uncertainty of 0.2 Hz. This uncer- 
tainty could best be lowered by illuminating the atoms 
with a larger beam and/or operating at lower values of 
the intensity. 

It is important to emphasize that the present uncer- 
tainties of these two shifts do not pose fundamental lim- 
its for future evaluations. Both the value of the shifts 
and the consequent errors scale linearly with lincwidths. 
In particular, an experimental linewidth of about 1 Hz 
could be produced with a reduced value of both the static 
magnetic field B = 0.3 mT and the probe light intensity 
I = 70 mW/cm 2 . In this case, with the present level 
of field calibrations the fractional frequency uncertainty 
would be less than 1 x 10~ 16 , and a feasible improvement 
by a factor of 10 in the calibrations would reduce the 
overall uncertainty to less than 1 x 10~ 17 . 

The frequency uncertainties due to lattice light have 
been studied in detail [ll|. The current systematic un- 
certainties (polarizability and hyperpolarizability) due 
to the optical lattice (magic wavelength X mag ic — 
759.3538 nm) at a depth of 500 E r are given in Table 
HI and considerable reduction of both is anticipated. 

TABLE I: Frequency uncertainty budget for the 174 Yb optical 
lattice clock. Some frequency shift values and their uncertain- 
ties depend on operating conditions and particular evaluation; 
typical values are given. 



Effect 


Shift (l(T io ) 


Uncertainty (10~ ib ) 


2 nd order Zeeman 


-18 - -36 


0.4 


Probe light 


6 - 12.0 


0.4 


Lattice Polarizability 


< 1 


0.6 


Hyperpolarizability 


0.33 


0.07 


Density 


-0.2 


1.0 


Blackbody shift 


-2.5 


0.25 


Syst. Total 


-15 --27 


< 1.5 



An important systematic effect to be accounted for is 
any density shift due to cold collisions between atoms in 
the ID optical lattice Q. This has been evaluated by 
reducing the total atom number captured in the 399 nm 
MOT and subsequently the lattice, changing the density 
in the range 0.5-1 xp where po ~ 10 11 cm -3 is the es- 
timated mean density in the lattice (see Fig. [T}:). The 
result, -0.1(0.6) Hz, is consistent with zero shift at this 
level of precision. The conservative uncertainty quoted 
in Table U is due to the small range of densities for which 
the system has been run reliably and the current opera- 
tion of the clock at the maximum atom number. Further 
reduction in the uncertainty of the density shift will oc- 
cur with better evaluations in a ID lattice or through 
the use of an under-filled 2D or 3D lattice, for which any 
such shift should vanish. 

Finally, the Stark shift induced by blackbody radia- 
tion (BBR) has been estimated by measuring the mean 
value of the temperature of the MOT chamber (295(3) K) 
and calculating the shift from the theoretical estimates 




Frequency (Hz) Frequency (Hz) 

FIG. 2: Frequency comparison of 174 Yb against a maser cal- 
ibrated to NIST-F1 (left) and Hg + (right) through use of 
the optical fequency comb 0, [l^] ■ Histograms and Allan de- 
viation generated from series of 1 s gate interval frequency 
counts. The linear fit on the two Allan deviations give a slope 
of 2xl0~ 13 r~ 1/2 and 5.5xl0" 15 r _1/2 respectively. The sta- 
tistical uncertainties for the two measurements are 0.9 Hz and 
0.10 Hz. 

in [12| . The uncertainty is limited by knowledge of Yb 
polarizability at the 10 % level, or 2.5 x 10 -16 fractionally. 
Precise measurements of the BBR shift or operation at 
cryogenic temperatures will be required to achieve sub- 
10 -17 uncertainty. 

Table U summarizes the systematic frequency shifts for 
magnetic field induced spectroscopy on 1 Sq ^-> 3 Pq transi- 
tion in 174 Yb. The values of the individual shifts depend 
on the particular operating conditions chosen (i.e., lat- 
tice power and frequency, bias magnetic field strength, 
and probe laser intensity). The quoted uncertainties of 
each frequency bias represent typical values. The actual 
uncertainties for a given comparison depend on the op- 
erating conditions and the proximity to, and quality of, 
systematic shift calibrations. With improvements in the 
experimental setup and diagnostics, reduced long-term 
variations in operation conditions will allow for less fre- 
quent systematic shift calibrations. 

Four absolute determinations of the clock transition 
frequency in 174 Yb were performed over several months 
against a maser calibrated with the NIST-F1 primary 
frequency standard [l3[ or against the optical Hg+ stan- 
dard using the optical frequency comb [3, [HI ■ As shown 
in Fig. [JU the noise contribution for both comparisons is 
mainly white in frequency, with a fractional instability 
of 2 x 10~ 13 t- 1/2 vs. the maser and 5.5 x 10~ 15 t" 1 / 2 
vs. Hg + . Recent improvements to the 578 nm stabiliza- 
tion cavity have resulted in an improvement of the Yb vs. 
Hg + fractional instability to less than 3 x 10~ 15 t~ 1//2 , 
with the stability reaching to below 1 x 10~ 16 in a couple 
of thousand seconds. For the comparisons against the 
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FIG. 3: Frequency evaluation of Yb clock transition. 
Two measurements were against Hg + (•) and the other 
two measurements were against a calibrated maser (a). 
The weighted mean of all the evaluations is z/( 174 Yb) = 
518 294 025 309 217.8(0.9) Hz (shaded area is the one a confi- 
dence interval). 

maser, statistical and calibration uncertainties are sig- 
nificant, and contribute a fractional uncertainty of about 
2.5 x 10~ 15 , depending on the length of the comparison. 
For comparisons against the optical Hg + standard, sta- 
tistical uncertainty can be made negligible (~ 2 x 10~ 16 ) 
with less than 20 minutes of averaging, and the current 
calibration of Hg+ to NfST-Fl is 0.65 x 10~ 15 . The ab- 
solute frequency for all the comparisons is then derived 
by correcting the statistical value with the calibration of 
the reference against NIST-F1 and for the gravitational 
shift caused by relative height differences. 

The final results of the frequency measurements of the 
clock transition in 174 Yb are reported in Fig. [3J Two 
values were obtained from optical comparisons against 
the Hg + standard and the other two were obtained 
though comparison against a N1ST-F1 calibrated maser. 
The weighted mean of all the evaluations is i/( 174 Yb) = 
518294025 309217.8(0.9) Hz. 

We have presented a frequency evaluation of the 1 Sq <-> 
3 Po clock transition in the 174 Yb isotope at a level ap- 
proaching that of the best atomic primary standards. 
Fractional frequency uncertainty below that of the NIST- 
Fl standard (~ 4 x 10~ 16 ) will be acheivable in the near- 
term with further measurements of the systematic fre- 
quency shifts. At this level, the current uncertainty in 
the room temperature BBR shift of Yb would become 
the limiting factor, and would need to be addressed for 
high accuracy ratio measurements against other optical 
standards [l], Q • Then straightforward modifications to 
the experiment (i.e. better vacuum to improve the life- 
time of the atoms in lattice, or multi-dimensional lattices) 
should lead to significant reductions in the other system- 
atic frequency shifts and uncertainties. Importantly, we 



foresee no significant barrier to achieving sub-10 _ accu- 
racy with an optical lattice clock based on magnetically 
induced spectroscopy of an even isotope of Yb. 
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